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Abstract 

The production of J/i|j pairs in proton-proton collisions at a centre-of-mass energy 
of 7 TeV has been observed using an integrated luminosity of 37.5 pb _1 collected 
with the LHCb detector. The production cross-section for pairs with both J/tp in 

the rapidity range 2 < < 4.5 and transverse momentum < 10 GeV/c is 

CT JAW =5.1 ±1.0 ±1.1 nb, 
where the first uncertainty is statistical and the second systematic. 
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1 Introduction 



The mechanism of heavy quarkonium production is a long-standing problem in QCD. 
An effective field theory, non-relativistic QCD (NRQCD), provides the foundation for 
much of the current theoretical work. According to NRQCD, the production of heavy 
quarkonium factorizes into two steps: a heavy quark-antiquark pair is first created pertur- 
batively at short distances and subsequently evolves non-perturbatively into quarkonium 
at long distances. The NRQCD calculations depend on the colour-singlet (CS) and colour- 
octet (CO) matrix elements, which account for the probability of a heavy quark-antiquark 
pair in a particular colour state to evolve into heavy quarkonium. 

Leading order (LO) calculations in the CS model [TJ El E] were first used to describe 
experimental data. However, they underestimate the observed cross-section for single 
J/ih production at high at the Tevatron [I] . To resolve this discrepancy the CO mech- 
anism was introduced [5]. The corresponding matrix elements were determined from the 
large-px data as the CO cross-section falls more slowly than the CS one. However, recent 
calculations [HI El E] close the gap between the CS predictions and the experimental 
data [10J reducing the need for large CO contributions. Thus, further experimental tests 
are needed. Pair production of quarkonium can cast light on this problem as this pro- 
cess depends heavily on the production mechanism. For both the CS and CO models, 
contributions from double parton scattering [HI [121 H3] could potentially be significant. 

The only observation of charmonia pair production in hadronic collisions to date was 
by the NA3 collaboration, who found evidence for J/ii) pair production in multi-muon 
events in pion-platinum interactions at pion momenta of 150 and 280 GeV/c [H] and in 
proton-platinum interactions at a proton momentum of 400 GeV/c [15]. The cross-section 
ratio a JAm/ a J/^ 

was measured to be (3 ± 1) x 10 4 for pion-induced production, where 
the inclusive J/ih production cross-section. At NA3 energies the main contribution 
to the J/4> pair cross-section arises from the quark-antiquark annihilation channel [TB] . 
This is not the case for proton-proton collisions at the LHC, where the gluon-gluon fusion 
process dominates p21 [18] . 

Theoretical calculations based on the LO production of CS-states predict that the 
total cross-section for J/ib-pair production in proton-proton interactions at y/s = 7 TeV 
is equal to 24 nb [121 ED]- These calculations take into account J/ibJ/ib, J/il)ih(2S') and 
il)(2S')ih(2S') production but do not include the possible contribution from double parton 
scattering. In the rapidity interval 2.0 < y J ^ < 4.5, relevant to the LHCb experiment, 
the expected value is 4 nb with an uncertainty of around 30%. At small invariant masses 
of the J/i[) pair a tetraquark state, built from four c-quarks, may exist [20] and would be 
visible as a narrow resonance in the mass spectrum. 

2 The LHCb detector and dataset 

The LHCb detector is a forward spectrometer [2TJ providing charged particle reconstruc- 
tion in the pseudorapidity range 1.9 < r] < 4.9. The detector elements are placed along 
the beam line of the LHC starting with the Vertex Locator, a silicon strip device that 
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surrounds the proton-proton interaction region. This reconstructs precisely the locations 
of interaction vertices, the locations of decays of long-lived hadrons and contributes to the 
measurement of track momenta. Other detectors used to measure track momenta com- 
prise a large area silicon strip detector located upstream of a dipole magnet with bending 
power around 4 Tm and a combination of silicon strip detectors and straw drift-tubes 
placed downstream. Two ring imaging Cherenkov detectors are used to identify charged 
hadrons. Further downstream, an electromagnetic calorimeter is used for photon and 
electron identification, followed by a hadron calorimeter and a muon system consisting 
of alternating layers of iron and chambers (MWPC and triple-GEM) that distinguishes 
muons from hadrons. The calorimeters and muon system provide the capability of first- 
level hardware triggering. 

The LHCb trigger system consists of three levels. The first level (LO) is designed 
to reduce the LHC bunch crossing rate of 40 MHz to a maximum of 1 MHz, at which 
the complete detector is read out. This is the input to the first stage of the software 
trigger, which performs a partial event reconstruction to confirm or discard the LO trigger 
decision. The second stage of the software trigger performs a full event reconstruction to 
further discriminate signal events from other pp collisions. To avoid that a few events 
with high occupancy dominate the CPU time, a set of global event cuts is applied on the 
hit multiplicities of each sub-detector used by the pattern recognition algorithms. These 
cuts were chosen to reject high-multiplicity events with a large number of pp interactions 
with minimal loss of luminosity. 

The data used for this analysis comprise an integrated luminosity of 37.5 pb _1 of pp col- 
lisions at a centre-of-mass energy of 7 TeV collected by the LHCb experiment between 
July and November 2010. This number includes the dead-time of trigger and data acqui- 
sition systems. During this period all detector components were fully operational and in a 
stable condition. The mean number of visible proton-proton collisions per bunch crossing 
was up to 2.5. 

The simulation samples used are based on the Pythia 6.4 generator [22] configured 
with the parameters detailed in Ref. [23]. The EvtGen [21] and Geant4 [25] packages 
are used to generate hadron decays and simulate interactions in the detector, respectively. 
Prompt charmonium production is generated in Pythia according to the leading order 
CS and CO mechanisms. 

3 Event selection and signal yield 

In this analysis the J/ih is reconstructed through its decay into a pair of muons. Events 
with at least four muons are selected. J/vJj> — > [i + \JT candidates are formed from pairs of 
oppositely-charged particles identified as muons that each have a transverse momentum 
greater than 650 MeV/c and that originate from a common vertex. Track quality is 
ensured by requiring that the Xt r / n df provided by the track fit is less than five. Well 
identified muons are selected by requiring that the difference in logarithms of the global 
likelihood of the muon hypothesis, provided by the particle identification detectors [26] . 
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with respect to the hadron hypothesis, A In £ M_h , be greater than zero. To suppress the 
contribution from duplicate particles created by the reconstruction procedure, if two muon 
candidates have a symmetrized Kullback-Leibler divergence (27] less than 5000, only the 
particle with the best track fit is considered. 

Selected M- + M-~ candidates with an invariant mass in the range 
3.0 < m^+n- < 3.2 GeV/c 2 are paired to form (\x + \i~)i(\i + \i~) 2 combinations. A fit of 
the four-muon candidate is performed [28J that requires the four tracks to be consistent 
with originating from a common vertex and that this vertex is compatible with one of 
the reconstructed pp collision vertices. To reject background where two J/4> candidates 
originate from different pp collisions, the reduced \ 2 °f this fit, xV n df, is required to be 
less than five. 

The number of events with two J/ih mesons is extracted from the single J/ih mass 
spectra. The invariant mass distributions of the first muon pair are obtained in bins 
of the invariant mass of the second pairj^] The single J/il> mass spectrum is modelled 
empirically using simulated events. This exhibits non-Gaussian tails on either side of the 
peak. The tail on the left-hand side is dominated by radiative effects in J/ip decay, while 
the right-hand side tail is due to non-Gaussian effects in the reconstruction. The shape of 
the distribution is described by a function that is similar to a Crystal Ball function [2H1 ED] , 
but with the power- law tails on both sides of the core Gaussian component. The position 
of the J/ih peak, the effective mass resolution and the tail parameters of this double- 
sided Crystal Ball function are fixed to the values determined from an analysis of the 
signal shape in the inclusive J/ih sample. Combinatorial background is modelled using an 
exponential function. This model is used to extract the yield of J/vp — > (\x + \xr) 1 in bins 
of the (\x + \x~)2 invariant mass. The extracted yield is shown in Fig. [T^i together with the 
result of a fit according to the model described above. The yield of events with double 
J/ih production given by the fit is N J ^ J ^ = 141 ± 19, where the statistical significance of 
this signal exceeds 6a. A fit with position and resolution of the signal peak left free was 
also performed and gave consistent results. 

Studies of single J/ih production indicate that the detector acceptance and efficiency 
is high for the fiducial range 2 < < 4.5 and < 10 GeV/c. The raw yield of 
events with both J/ih mesons within this range is 139 ± 18. The yield of events with both 
J/ih mesons in the fiducial range and explicitly triggered by one of the J/vp candidates 
through the single muon or dimuon trigger lines [31], is found to be 116 ± 16. This sample 
is considered for the determination of the production cross-section. 

The contribution to the yield from the pileup of two interactions each producing a sin- 
gle J/ih meson is estimated using simulation together with the measured J/ih production 
cross-section [32]. This study shows that for the 2010 data-taking conditions the back- 
ground from this source is expected to be less than 1.5 events and hence can be neglected. 
In a similar way the contribution to the yield from events with J/vp mesons originating 
from the decays of beauty hadrons is found to be negligible. 

2 The M. + M-~ pair with lower transverse momentum is chosen to be the first pair. 
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"W-) 2 [GeV/c 2 ] m (v+v-) 2 [GeV/c 2 ] 

Figure 1: The fitted yields of J/vJj — > ([i + \x~) 1 in bins of ( (x — ) 2 invariant mass: (a) the 
raw signal yield observed in the data; (b) the efficiency corrected yield (Sect. [6]). The re- 
sult of a fit with a double-sided Crystal Ball function for the signal and an exponential 
background is superimposed. 



4 Efficiency evaluation 

The per-event efficiency for a J/ib-pair event, ^j/^j/^, is decomposed into three factors, 

tot _ reco AtlD trg /,\ 

£ J/*J/4> _ £ JAWA1> x £ J/t|)J/t|) x e JAP JAP ' \ L ) 

where Sja^jaj, is the product of the (geometrical) acceptance with reconstruction and 

selection efficiency, is the efficiency for muon identification and ^j/^j/^ is the trigger 

efficiency for reconstructed and selected events. 

The efficiency for the acceptance, reconstruction and selection for the two J/i[> mesons 
is factorized into the product of efficiencies for the first and second J/ih , 

reco ._ reco / J/^i „.J/4» i I pntJ ,Q* M v c- rcco (n^ 2 ti J /^2 I nne i9* A 

£ JA|)J/t|) — £ JAl> '2/ , I cosvjy^JI x yp T ,y 2 , | cos^j/q,^ I . I^zj 

The single JAJj efficiency e^ is a function of the transverse momentum £>t ; rapidity y 
and | cos i?* |, where $* is the angle between the \x + momentum in the J/4> centre-of-mass 
frame and the J/vp flight direction in the laboratory frame. It is evaluated using simulation. 
The validity of the factorization hypothesis of Eq. ^ is checked with simulation and based 
on these studies a correction factor of 0.975 is applied to Sj^jaj, • For the simulated data of 
single prompt J/tp production the cut on the muon likelihood is not applied and that on 
X 2 (JAP J/ih)/ n df is replaced with a similar cut on the single J/vp , x 2 (J/4')/ n df < 5. The 
efficiency of the cut on x 2 /ndf is estimated from the data and compared to the simulation. 
Based on these studies a correction factor of 1.026 is applied to £j/^°t/^, and a systematic 
uncertainty of 3% is assigned. The efficiency is also corrected by a factor 1.024 ± 
0.011 [32], that accounts for the ratio of the reconstruction efficiency of the muon detector 
observed in the data compared to the simulation. 
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The muon identification efficiency is extracted from the analysis of the inclusive 
J/ih sample. Two efficiencies are evaluated: the single muon identification efficiency 
e^ ID and the J/4> efficiency e^*. Both are measured as a function of the value of the 

cut made on Aln£ M ~ h . The squared efficiency (e^ ) 2 and are found to be equal to 

better than one per mille. The value of = ( e j/™) 2 = (91.0 ± 0.1)% has been used 

as a global factor for the evaluation of the total efficiency using Eq. 0. As a cross-check, 
the efficiency of the muon identification for J/4> pairs has been estimated from the signal 
itself. Though statistically limited, the value is consistent with that given above. 
The trigger efficiency is calculated to be 

e^u =l-(l-e t * (p^\v^A) x (l-e^. (vf>\v^ 



j/* j/* = 1 - 1/ - £ J/4> [Pt ,V'^) ) x " £ jfi> [Pr ',V 

The trigger efficiency for a single J/ih , , is determined directly on data from the inclu- 
sive prompt J/ih sample as a function of pp and rapidity y. The efficiency is determined by 
classifying an event which would also have been triggered without the J/ij> as TIS (Trigger 
Independent of Signal), and/or classifying the event where the J/4> alone is sufficient to 
trigger the event as a TOS (Trigger On Signal) event [33l 34J. The LHCb trigger system 
records all the information needed for such classification. Events can be classified as TIS 
and TOS simultaneously (TIS&TOS), which allows the extraction of the trigger efficiency 
relative to the off-line reconstructed and selected events from data alone 

trg _ iV TISfcTOS 
£ -VA> ~ jyTis ' 

where N Tls is the number of TIS events, and _/v TIS&:TOS is the number of events that are 
simultaneously TIS and TOS. The method has been cross checked using Monte Carlo 
simulation. 

The effect of the global event cuts applied in the trigger has been studied in detail 
for inclusive J/ih events in Ref. [32]. Since the subdetector hit multiplicity observed in 
single and double J/4> events is similar, the efficiency of the global event cuts, (93 ±2)%, 
is taken from that study and applied to ^j/^j/^- 

For selected J/ib-pair events the mean value of ^j/^j/^ is 31% and it varies from 0.8% 
to 70%. The mean value for e^,^ is 85% and it varies from 61% to 93%. 



5 Systematic uncertainties 

Systematic uncertainties affecting the cross-section measurement have been evaluated 
properly taking correlations into account where appropriate. The dominant source of sys- 
tematic uncertainty is due to the knowledge of the track-finding efficiency. An uncertainty 
of 4% per track is assigned based on studies comparing the reconstruction efficiency in 
data and simulation using a tag and probe approach [34J. 

The second major source of uncertainty is due to the evaluation of the trigger efficiency. 
The method discussed in Sect. [4] has been cross-checked in several ways, in particular, by 
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using events triggered by the first or second J/ip only. The observed differences lead to 
the assignment of an 8% systematic uncertainty. 

A further source of uncertainty is the determination of the per-event efficiency defined 
by Eq. Q. This is estimated to be 3% by varying the uncertainties of the various factors 
entering into Eq. 

The unknown polarization of J/ih mesons affects the acceptance, reconstruction and 
selection efficiency £j^° [22] • In this analysis the effect is reduced by explicitly taking 
into account the dependence of the acceptance on e T f^ on |cos i?* | in the efficiency de- 
termination (Eq. [2]). The remaining dependence results in a systematic uncertainty of 
5% per J/ii) . 

Additional systematic uncertainties arise due to the difference observed between the 
data and simulation for the behaviour of the cut on x 2 (3%), the global event cuts (2%), 
and uncertainty of 1.1% per J/ih associated with the efficiency for muon identification, 
as discussed in Sect. [4} The systematic uncertainties associated with the other selection 
criteria and the J/ih lineshape parametrization are negligible. 



Table 1: Relative systematic uncertainties on the cross-section measurement. The total 
uncertainty is calculated as the quadratic sum of the individual components. 

Source Systematic uncertainty [%] 

Track finding efficiency 4x4 

Trigger efficiency 8 

Per-event efficiency 3 

J/vp polarization 2x5 

Data/simulation difference for % 2 /ndf 3 

Global event cuts 2 

Muon identification 2 x 1.1 

Luminosity 3.5 

J/iL) — > (J. branching ratio 2x1 

Total 21 



The luminosity was measured at specific periods during the data taking using both 
van der Meer scans [35] and a beam-gas imaging method [36]. The instantaneous lumi- 
nosity determination is then based on a continuous recording of the multiplicity of tracks 
reconstructed in Vertex Locator, which has been normalized to the absolute luminosity 
scale. Consistent results are found for the absolute luminosity scale with a precision of 
3.5%, dominated by the beam current uncertainty [3T| I3"8"]. 

The relative systematic uncertainties are summarized in Table [TJ where the total 
uncertainty is defined as the quadratic sum of the individual components. 
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6 Cross-section determination 



The efficiency-corrected yield for events with both J/xJj candidates in the fiducial region is 
extracted using the procedure discussed in Sect. [3] To account for the efficiency a weight 
u, defined as 

, ,-1 _ ptot 

where ^j/^j/^ is the total efficiency defined in Eq. (jl]), is applied to each candidate in 
the sample. 

The corrected yield of J/ih — » in bins of ( ) 2 invariant mass is shown in 

Fig. [TJo. As previously described, to extract the yield a fit with a double-sided Crystal 
Ball function for the signal, together with an exponential function for the background 
component, is performed. Again, the position of the J/ih peak and the effective mass 
resolution are fixed to the values found in the inclusive J/ip sample. The event yield after 
the efficiency correction is 

N#W = 672 ± 129. 

The cross-section for double J/ih production in the fiducial range 2 < y J ^ < 4.5 and 
< 10 GeV/c is computed as 



rcorr 



JAW _ JAW 



where C = 37.5 ± 1.3 pb 1 is the integrated luminosity and £y+|x- = (5.93 ± 0.06)% [39 
is the J/ih — > branching ratio. The result is 

CT JAW =5.1 ±1.0 ±1.1 nb, 

where the first uncertainty is statistical and the second systematic. 

Using the measured prompt J/ih production cross-section in the same fiducial re- 
gion |32j and taking into account the correlated uncertainties, the ratio of cross-sections 
oWW/aW i s calculated to be 

oWW/eW = (5.1 ±1.0 ±0.6^) x 10" 4 , 

where the first error is statistical, the second systematic and the third is due to the 
unknown polarization of the prompt J/ih and J/ih from pair production. 

The differential production cross-section of J/ip pairs as a function of the invariant mass 
of the J/ih J/ih system is shown in Fig. |2j The whole analysis chain has been repeated 
for each bin of J/ipJ/ih invariant mass to get the differential production cross-section. 
The bulk of the distribution is concentrated in the low invariant mass region. A theoretical 
prediction for the shape of this distribution taking into account both direct production and 
feeddown from ty(2S) decays [20] is overlaid. Within the available statistics the agreement 
between the data and the prediction is reasonable. 
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mj/^j/^ [GeV/c 2 ] 

Figure 2: Differential production cross-section for J/i|) pairs as a function of the invariant 
mass of the J/i|)J/i|) system. The points correspond to the data. Only statistical uncer- 
tainties are included in the error bars. The shaded area corresponds to prediction by the 
model described in Ref. [20] . 

7 Conclusions 

The production of J/ih pairs in proton-proton collisions at a centre-of-mass energy of 7 TeV 
has been observed with a statistical significance in excess of 6a. The data are consistent 
with the predictions given in Refs. [191 EQ]. The higher statistics that will be collected 
during the 2011 data-taking period will allow the kinematic properties of these events to 
be studied and different production models to be probed. 
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